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Abstract — Due to hardware and computational constraints, 
wireless sensor networks (WSNs) normally do not take measure- 
ments of time-of-arrival or time-difference-of-arrival for range- 
based localization. Instead, WSNs in some applications use range- 
free localization for simple but less accurate determination of 
sensor positions. A well-known algorithm for this purpose is the 
centroid algorithm. This paper presents a range-free localization 
technique based on the radical line of intersecting circles. This 
technique provides greater accuracy than the centroid algorithm, 
at the expense of a slight increase in computational load. 
Simulation results show that for the scenarios studied, the radical 
line method can give an approximately 2 to 30% increase in 
accuracy over the centroid algorithm, depending on whether or 
not the anchors have identical ranges, and on the value of DOI. 

Index Terms — Wireless sensor networks, radical line, localiza- 
tion algorithm, centroid algorithm. 

I. Introduction 

A wireless sensor network (WSN) [T) typically consists 
of anchors and sensors communicating with each other. An 
anchor broadcasts its position coordinates, together with oper- 
ating instructions, to the sensors. A sensor needs to determine 
its position to report to the anchors. Position determination can 
come from time-of-arrival, time-difference-of-arrival or angle- 
of-arrival measurements [2|. But when the sensors are low 
cost, low power and expandable units, with limited resources 
for computation, they often rely on range-free (RF) localiza- 
tion instead 0. 

In RF localization, a sensor P determines its unknown 
position P = [x,y] T from N in-contact anchors <Zj at known 
positions a.; = [xi : yi] T and having radio ranges R%,i = 
1, N. The sensor position must satisfy 

\\P-Oif' 2 = [{x-x i )* + { y -y i f\W 
^Ri, i = l,2,...,N. 

Solving (JTJ requires nonlinear programming, and there is not 
a unique answer. The Centroid Algorithm (CA) [3] gives a 
simple estimate P = [x, y] T , where 
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But P from (fj) sometimes is outside the region of intersections 
(RI) of the circles centered at a, with radii R, as defined by 
dTJ. For example, the P in Fig. 1 is outside the RI (shaded 
area) of the three circles. 




Fig. 1. The CA and radical line solutions 



This paper proposes a new RF algorithm that has better 
accuracy than CA, but with a marginal increase in computa- 
tions. However, the additional computations are well within 
the capability of present day sensors. 

The line joining the intersection points of two circles is the 
radical line (RL) [4|. In Fig. 1, the RI contains a segment of 
the RL of any two circles, and the three RLs meet at a point 
inside the RI. Indeed, |4] proves that for three circles whose 
centers are not collinear, their three RLs always intersect at a 
point. Although sometimes this point can be outside the RI, it 
is inside in most cases. 

In the following, Section II gives the development of the 
RL algorithm (RLA). Section III contains simulation results, 
which show that the RLA is more accurate than the CA, 
especially when the radii Ri are different. Conclusions are 
given in Section IV. 

II. The Radical Line Algorithm 

In WSNs, a sensor can determine whether it is in the 
transmission range of an anchor node according to the beacon 
signal received from the anchor. Most literature on range-free 
localization assume a nominal range (or detection range) R, 
i.e., an anchor can communicate with a sensor within R meters 
from it. However, the actual range in practice is dependent on 
the propagation conditions. A measure of the variation in range 



coverage is the degree of irregularity (DOI). Its value denotes 
the maximum range variation per unit degree change in the 
direction of radio propagation. Recently, [5| gives a condition 
required for a successful anchor-to-sensor contact. Let W(oj) 
be the power received by a sensor from a%, Q be the ambient 
noise power, and S the interference power in the WSNs. Then 
there is a contact only if 

W(a., 



Q + S 



> TH, 



(3) 



where TH is a hardware dependent threshold. 

Let a sensor P be at an unknown position P = [x, y] T , in 
contact with N anchors ai at known positions = [xi,yi] T 
and having radio ranges Ri. Hence P must lie in the RI of 
the N circles, centered at ai with radii Ri. Depending on N, 
there are three cases to consider. 



A. N > 3 

For N circles, there are 21(^-2)! RL S - To reduce compu- 
tations, RLA selects only the RLs of the two circles whose 
centers are separated by the largest distance among the N 
circles. The idea behind this choice is that the RL of these 
two circles will be the shortest, and hence their RL has the 
highest probability of appearing inside the RI of all the N 
circles. 

Let 

„ 11 1/2 



= dji, i,j = l,...,N 



(4) 



be the distance between the centers of eij and aj and let d q t be 
the maximum of the values in ©. For illustration simplicity, 
let q = l, and k = 2. Referring to Fig. 2, the end points of the 
RL are I a = [x a ,y a ] T and I b = [x b ,y b ] T , and O = [x ,y ] T 
is the intersection between the RL and the line joining a% and 
a 2 . 
Let 



and 



doi = \\0 - ai 
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It follows that 
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and 



di 2 +m' = R4. 
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Subtracting dH) from (O gives 

2(x 2 - xt)x + 2(y 2 - yi)y a =Rf-Rl + k 2 - h, (9) 
where 



ki = x t + y t 



(10) 



Further, equating the slopes of the a\ to O and a 2 to a\ lines 
in Fig. 2 yields 

2/2 - Vo 2/2 - 2/1 
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X2 — X\ 



(11) 



giving 

(2/2 - yi)x Q - (x 2 - xi)y D = x 2 (y 2 ~ 2/1) - 2/2(^2 - %%)■ (12) 



Solving (O and (fT2l then gives O(x ,y ). 
Let di2 = D. Then 



R\-d 2 ol =m 2 = R%-{D-d xf, 



so that 
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and 



m = (R\-d 2 ol ) 1 ' 2 . (15) 
Now in Fig. 2, the following trigonometric relations hold: 

x -x a 2/0-2/1 

= — (16) 



and 
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From ( TToT l and dT7l >. the coordinates for I a are 



Xa = X - -T—{Vo - 2/l) 
«ol 



and 
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Va = 2/o + -J—(Xo ~ X\). 
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Following the same procedure gives 
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x b = x + —{y - 2/1) 
dot 



Vb = y - -i—(x ~ x x ). 
dot 
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Next, RLA selects L test points t\ = [xi,yi] T , I = 
1, . . . , L, on RL, by taking equal increments between I a and 
l h to give 

l(x b - x a ) 



( ? ) and 
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L is a user parameter, depending on the resolution required. 
In the simulation experiment in Section EH, L — 4. At each 
ti, RLA checks whether ti is inside the RI, and if not, how 
far away from the RI it is, by computing the error 
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and then summing the errors over all a, to give 
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Fig. 2. The end points of an RL 



If an Si = 0, the corresponding t\ is inside the RI and is the 
estimate for P. If all Si > 0, the RL is not inside the RI of 
the N circles. It is then necessary to compute the CA errors 
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where e C i comes from d24l i. with c = [x, y] T from (ffji replacing 
ti. The final estimate for P, P, comes from choosing the ti 
or c, whose corresponding Si or S c is the minimum. 

B. N = 2 and N = 3 

When N = 2, P is the same as 0{x o ,y o ). When N = 
3, RLA computes the intersection of the three RLs. Let that 



intersection point be I 
circles yields 



[xi, yi] . Extending Fig. 2 to three 



(xi - Xi) 2 + ( yi - y t ) 2 + h 2 
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where h 2 < m 2 . Subtracting this expression for i = 2, 3 from 
that for i = l results in 
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Solving d28T i gives 
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If the determinant of A equals 0, then the three circles are 
collinear. Or if ||7 — a.^ 1 1 > R% for any i, then / is outside 
the RI. For these two cases, RLA takes the centroid of the two 
circles with the largest separation as P. 



III. SIMULATION RESULTS 

In the simulation experiments, the WSN has an area of 100 
m x 100 m, and contains 100 randomly placed (different for 
each trial) sensors. For a given number of anchors (NA), occu- 
pying random (different for each trial) but known positions, the 
number of anchors in contact with an arbitrary sensor can 
vary from 2 to NA. Some anchors have R = R max — 45m, 
and some have R = 0.5R max . The localization errors decrease 
with increasing NA. For 100 independent trials, the error as a 
fraction of R m ax is 
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In (|32l , Pj(i) is the true z'th sensor position at trial j, and pj(i) 
is its estimate. 

In an experiment where DOI = 0, a sensor that lies 
within the nominal i?j of an anchor is in contact with that 
anchor. When DOI ^ 0, the actual Ri is smaller, given by 
Ri(DOI) = Ri(l - DOI). 

Fig. 3 plots e(NA) for both CA and RLA with all anchors' 
R = R m ax and as NA varies from 24 to 36, at a DOI = 0. 
The results show that RLA has lower e(NA) than CA. 

Fig. 4 plots e(NA) for both CA and RLA and all R = R max 
as NA varies from 24 to 36, at a DOI = 0.1. The results show 
that RLA has lower e(NA) than CA. 

Fig. 5 plots e(NA) for both CA and RLA with different 
transmission ranges, i.e., some cij have R — R m ax — 45m 
and some have R — 0.5R max as NA varies from 24 to 36, 
at a DOI = 0. The improvement of RLA over CA is more 
significant than when all anchors have R — R max . 

Fig. 6 plots e(NA) for both CA and RLA with different 
transmission ranges as NA varies from 24 to 36, at a DOI = 
0.2. A comparison of the errors in Figs. 4-6 reveals that 
RLA has increasing accuracy over CA, when DOI increases. 
The improvement is more significant when the anchors have 
different ranges. 

Fig. 7 is a snapshot of one trial in the anchor-sensor 
geometry with NA = 30 and different transmission ranges, 
together with the placement of P. A dotted line joins P to P. 
Comparing Fig. 7(a) to Fig. 7(b), the dotted lines for RLA are 
generally shorter than those for CA. 

Fig. 8 plots e(NA) for both CA and RLA as the number 
of sensors varies from 50 to 80 with the same transmission 
range, at DOI = 0.1. The number of anchors NA = 30. 

Fig. 9 plots e(NA) for both CA and RLA as the number 
of sensors varies from 50 to 80 with different transmission 
ranges, at DOI = 0.1. The number of anchors AM = 30. 
Comparing Fig. 8 with Fig. 9 shows that the accuracy gain 
of RLA over CA is higher when the anchors have different 
ranges, than when they have the same range. 

IV. Conclusions 

Range-free localization, while not as accurate as range- 
based, has the principal advantage of simplicity, i.e., there is 
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Fig. 4. The average localization error vs. the number of anchors (DOI=0. 1 
and the same transmission range) 



Fig. 6. The average localization error vs. the number of anchors (DOI=0.2 
and different transmission ranges) 



no requirement for special hardware to measure time-of-arrival 
or time-difference-of-arrival. This is important for WSN in 
which sensors are low cost units, and in some applications 
where knowing accurate sensor positions is not critical. While 
determining P from CA is simple, it is possible to improve 
on its accuracy with some additional computations. The RLA 
provides such an option and the simulation results in Section 
III show that there is approximately a 2 to 30% gain in 
accuracy, depending on whether or not the anchors have 
identical ranges, and on the value of DOI. 
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Fig. 8. The average localization error vs. the number of sensors (DOI=0. 1 
and the same transmission range) 
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(a) Localization error of RLA (DOI=0.1, error = 0.1929) 
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(b) Localization error of CA (DOI=0.1, error = 0.2872) 
Fig. 7. Location error with different transmission ranges 
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Fig. 9. The average localization error vs. the number of sensors (DOI=0. 1 
and different transmission ranges) 



